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The interplay between flavour symmetries connecting B s — > KK decays 
' with the recently measured B4 — > K°K° decay and QCD Factorisation 

opens new strategies to describe the decays B s — > K°K° and £? s — * K + K~ 
in the SM and in supersymmetry. A new relation, emerging from the sum- 
rule for the -B s — *■ K°K° decay mode, is presented offering a new way to 
determine the weak mixing angle <ps of the B s system. 

o 

PACS numbers: 13.25.Hw,11.30.Er,11.30.Hv,12.39.St 

1 , 1. Introduction 

The huge effort on the experimental side at present B facilities (Babar, 
Belle and CDF) to increase the precision on data measurements force us to 
revise the strategies on the theory side to produce more accurate predictions. 
Non leptonic B decays offer different strategies to determine the Unitarity 
Triangle, to search for New Physics (NP)[T] but also to rule out models [21 • 
While a lot of attention has been devoted to the B — > 7rif[31[HG3[Hl[7] decay 
modes, here we will focus on B —* KK decays that has been observed at 
CDF |5] (B s -> K+K~) and at Babar [9] and Belle [TO] (B d K°K°). 

There are two main approaches in the literature to describe B — > KK 
decays: flavour symmetries and l/m&— expansion methods (QCD Factori- 
sation [Til 12]/ soft collinear effective theories [13] or PQCD [13]). Each 
of those methods has pros and cons, that we will discuss in turn. Flavour 
symmetries, like U-spin symmetry that relates B s — ► K + K~ with 2?^ — > 
7r + 7r _ [151 (El El [E]) provide a model independent analysis and extract 
most of the needed hadronic parameters from data. However, this method 
has the disadvantage that it relies strongly on the accuracy of data, and, at 
present, there is still not full agreement between Babar and Belle data on 
the CP asymmetries of the B^ — * tt + tt~ mode. As a consequence, error bars 
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are still quite large, see for instance, the prediction BR(S S — * K + K ) = 
( 35 ±2o) x 10_6 [7\ or > more recently, 4.2 x KT 6 < BR(B S -> K+R-) < 
60.7 x 10~ 6 [18]. Also when relating B d -> 7r+7r~ with B s -> K + K~ some 
of the needed U-spin parameters can only be roughly estimated and they 
are usually taken to be of the order of 20%. 

Concerning l/m^-expansion methods, here we will focus on QCD Fac- 
torisation (QCDF) [111 112} [19]. The main idea is to exploit the existence of 
a large scale m& S> Aqcd together with colour transparency, that applies 
when the outgoing meson that does not contain the spectator quark is very 
energetic. At leading power in A/ m& ah long distance contributions can 
be parametrized in terms of form factors and light cone distribution ampli- 
tudes, while the contribution from energetic gluons comes in a perturbative 
series in a s and it is incorporated into the hard scattering kernels. QCDF 
predicts some of the hadronic parameters reducing the error bars, however 
in the computation one has to face chirally enhanced IR divergences. They 
are formally suppressed by a power of 1/rrib, but can be numerically sig- 
nificant. They are modelled and induce an important uncertainty in the 
predictions. 

However, there is a third possibility and it is the proposal presented 
in [20] that combines QCDF and Flavour symmetries giving rise to rather 
accurate predictions for the branching ratios of the above mentioned B s — > 
KK decays in SM and in super symmetry. Moreover, the method predicts 
some of the SU(3) breaking parameters which can be useful for other flavour 
approaches and, at the same time, deals with the problem of the chirally 
enhanced IR divergences coming at order A/m^ (see also [21]). 

Since IR divergences play a central role in this discussion, it is worth to 
mention the two sources of IR divergences in QCDF: 

• Hard spectator-scattering: Hard gluons exchange between spectator 
quark and the outgoing energetic meson gives rise to integrals of the 
following type (see [12] for definitions): 



Hi(MiM 2 ) = C Cdx f\y 
Jo Jo 



xy xy 
where the second term (formally of order A /mj) diverges when y — * 1. 

Weak annihilation: These type of diagrams also exhibit endpoint IR 
divergences as it is explicit in the corresponding integrals: 

A\ = na s [ dxdy l<f>M 2 (x)®M 1 (y) n 1 

Jo I Vy(l-xy) 

+ 4-1 + r x lr x 2 • 
x A y\ xy ) 
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Both divergences are modelled in the same way 
1 ^* mi (y) = * mi (l) /* 5+ t ^[$ mi (y)-$ mi (l)l =<& mi (l)JT^+/- 



o y ./o y 



where r is a finite piece and the divergent piece is cut-off by a physical scale 
of order Aqcd with an arbitrary complex coefficient to take into account 
possible multiple soft scattering: Xh A = (1 + Ph,a) ln(mfi/A). 



2. Sum-rules: a and </> s 

The SM amplitude for a -B-decay into two mesons can be split into tree 
and penguin contributions :A = A{B q MM) = Ai 9) T|f + \ { c q) Pjf , 
with C denoting the charge of the decay products, and the products of CKM 
factors Ap 9 ^ = VpbV*„. Tree and penguin contributions in B d — > i^°-ftT° in 
QCDF are: 

T'dO W ^ „ « | fl« I r) 1 OH /Oil 

— a A — -a iEW + p 3 + — -PsEW ~ Paew 

P d0 = a% — -a% EW + 13% + 2fil — -fil EW — PIew i (1) 

where P d0 = P d0 /A d KK , f d0 = T d0 /A d KK , the super-scripts identify the 

channel, the normalisation is A q KK = M% Ftf q ~* K (p)f K G F /y/2 (see [121120] 
for the corresponding expression of the B s — > KiT channels). Following 
the observation in [20] that the structure of the IR divergences is the same, 
independently of the charm or up quark running in the loop, we identified 
an IR-safe quantity at NLO in QCDF that we called A d = T d0 - P d0 . 
All chirally enhanced IR divergences cancel exactly in this quantity at this 
order. Its explicit expression in terms of the coefficients in Eq.([T]) is: 

A d = A d kk [at -a% + % - (3 C 3 + 2$ - 2/3 4 c ], 

where electroweak contributions are neglected. This quantity can be safely 
evaluated in QCDF and the result found in [20J was: A d = (1.09 ± 0.43) • 
10~ 7 + i(— 3.02±0.97) • 10~ 7 GeV. The largest uncertainty entering comes 
from the ratio m c /m& and the scale dependence. Interestingly, this quantity 
can be expressed in terms of observables, providing a relation between the 
direct induced CP-asymmetry (A~9 ), the mixing induced CP-asymmetry 
(Af^ ix ) and the branching ratio (BR d0 ) of B d -» K°K°: 

BR d0 

|A d | 2 = — — {xi + [x 2 sm(f) d - x 3 cos(f) d ]A™ ix - [x 2 cos (p d + x 3 sin^]^ } , 



1 Conventionally, we will call "tree" the piece proportional to A^ 9 ' and "penguin" the 
piece proportional to \i q ^ , even if applied to decays with no actual tree diagram. 
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where \A^\ 2 + l^irl 2 + l^mixl 2 = 1> 4>d is the weak mixing angle for the 
Bd system, = r^y Mg d — 4M|- / (327rMg d ) and ccj are functions of 7 and 
CKM elements. All SM inputs are taken as in [20J following |22j. 

Moreover, it was found in [23] that this sum-rule encodes also a very 
interesting information. It provides a new way of measuring sin a: 



(See |24| for a recent review on the extraction of a). In a similar way, it was 
found in [20] a corresponding IR safe quantity A s = T s0 — P s0 and sum-rule 
for the decay B s K°K°: 



A s | 2 = — — {yi + [y 2 sin 4> s - y 3 cos 4> s ]A s £ ix - [y 2 cos 4> s + y 3 sin 4> S ]A S £} , 



which provides a completely new way to determine the weak mixing angle 
<j) s that we present here (see [20J for definitions): 



This implies that a measurement of the branching ratio, direct CP asym- 
metry and mixing induced CP asymmetry of the decay B s — ► K°K° auto- 
matically translates into a value for sin 2 4> s /2. Finally, given the relation 
A s = /Ad, where / = A S KK / ' A d KK a new relation between sin a and sin<^ s /2 
immediately emerges. 

3. Description of the Method: Flavour Symmetries & QCDF 

The SM amplitude of this b — > s penguin decay is given by: 



Its dynamics is described in terms of three parameters: \T \, \P s0 \ and 
the relative strong phase arg(i- >s0 /T s0 ) (remember that T s0 stands for the 

(s) 

piece proportional to A« but it is not due to an actual tree diagram in this 
case). Its hadronic parameters can be related via U-spin with the hadronic 
parameters (|T |, \P d0 \ and aig(P d0 /T d0 )) of the also penguin governed 
mode Bd — > K°K°. This has several advantages: first, we can expect 
similar final state interactions (although not equal), second, the sources of 
U-spin breaking can be better controlled using QCDF. These sources are: 





A = A(B S -> K°K°) = Al s) T s0 + X^P 
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i) the factorisable ratio / = A s KK j A d KK (extrapolated from the lattice) ii) 
U-spin breaking l/m& suppressed terms 5cti and Sfy: sensitive to the differ- 
ence of B d and B s distribution amplitudes and spectator quark dependent 
contributions coming from a gluon emitted from the d or s quark. This 
leads to the relations \P s0 /(fP d0 ) - 1| < 3% and \T s0 /(fT d0 ) - 1| < 3%. 

The next step is to determine the hadronic parameters (T d0 , p d0 ^ of the 
decay B d — * K°K°. This is done using as inputs the BR(B d — ► K°K°) and 
A d (from QCDF). The direct CP asymmetry A dir (B d K°K°) (denoted by 
^cfcr) WU ^ ^ e taken as a free parameter. The combination of those constraints 
gives rise to a set of non-linear equations (see definitions in |20j): 

x c + iyc = -Ad(l - cos j /R) /a, 

r 2 = pg/[a|AW| 2 ]-[sin 7 |A d |/(ai?)] 2 , 
?/pxa = ?/axp - ^ r /(2|A(f)AW| sin 7 ) , (2) 

that determines P rf0 = xp + iyp, then using one gets T rf0 . Two remarks 
are important here: first, there is a twofold ambiguity in the sign of ImP d ® 
(solved in the next subsection). Second, current data together with our 
knowledge on A d limits the Af^ r asymmetry (by means of Eqs.(|2])) within 
a restricted range between -0.2 < Af ir < 0.2. 

Finally, our SM predictions for the branching ratio and CP asymmetries 
of B s — ► K°K° are obtained using the relations between P s0 *-> P d0 and 
T s0 «-> T d0 mentioned above and including all U-spin breaking sources to- 
gether with the QCDF uncertainties in A d . The resulting predictions are 
BR(B S -> K°K°) = (18 ± 7± 4db (2)) x 10~ 6 , \A dir (B s -> i^ ^ )] < 
1.1 x 10~ 2 and | Anises -> Er K°)| < 1.5 x 10~ 2 . The sign of these asym- 
metries can be fixed once A^ r will be measured with enough accuracy. 

3.2. B s K+K- 

The analysis of the decay mode B s — > K + K~ follows similar steps, but 
with some important differences: i) we will use U-spin and isospin to connect 
B d -> K°K° with B s -> K + K~, ii) £ s -> K + K~_ contains a tree (denoted 
by ai in QCDF) with no counterpart in S d -» if °iT , iii) Si?(5 s -» K+R-) 
has been measured with excellent precision at CDF [8] and iv) we will use 
the information only on the sign predictions (not the absolute value) for 
the CP-asymmetries of B s — > K + K~ from the strategy that uses U-spin to 
relate B s — ► K + K~ with — > 7t + 7t~. 

The hadronic parameters describing the amplitude: 

A(B S -» K + K~) = \$T S± + A( s) P s± 

are obtained from the relations containing the sub-leading l/m& U-spin 
breaking: |P s± /(/P d0 ) - 1| < 2%, |r a± /(A£ fc ai) - 1 - T d0 /(^ fc ai)l < 4%. 
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Those errors, estimated within QCDF, are stretched roughly by a factor 
two to be conservative. 

The two-fold ambiguity on the sign of ImP d0 is lifted here using U-spin 
arguments based on the B^ — ► 7r + 7r~ strategy. While the signs of ImP and 
A s d - are correlated (being both positive or negative), the prediction based 
on the Bd — » ir + ir~ strategy points towards a positive sign for A s j^_ |18j . 
discarding then the solution with ImP d0 < 0. 

Our result in the SM for the branching ratio of B s — > K + K~ averaging 
over all values of Af r is BR(B S -> K+K~) = (20 ± 8 ± 4 ± (2)) x 10~ 6 , 
where the last error in parenthesis stands for a rough estimate of finite, 
non-enhanced A/m& corrections. Finally, confronting our predictions for 
B s -» K+K- with the data on B d -> tt+tt' [25]: |T^±| = (5.48±0.42) x 10~ 6 

and P^I T H = °- 13 ± °- 05 > arg(p^/T^) = (131 ± 18)°, provides a 
double information. First, we can give predictions for the U-spin breaking 
parameters: K c = \T S± /T d ±\ = 2.0 ±0.6 and £ = \P S± /T s± \/\P d ± /T d ±\ = 
0.8 ± 0.3 connecting S s — > K + K~ with 1?^ — ► 7r + 7r~. These parameters 
can be compared with the QCD sum rules predictions in ref. Notice 
that while QCD sum rules gives only the factorizable part, our predictions 
include, in principle, the full contribution. Second, a comparison between 
the two relative strong phases arg (P s± /T s± ) and arg \P d ^ /T d ^j selects 

A^ r > 0. Then, if we restricts only to positive values of A^ according to 
the previous arguments, our SM predictions turn out to be [20] |2"7]: 

BR(B S ^K + K~) = (17±6±3± (2)) x 10~ 6 , 
-0.22 < A s d f r < 0.49 and -0.55 < A s ± x < 0.02. (3) 

Another argument in favour of A$ > comes from the preference of A^ < 
of the U-spin based B^ — ► tt + tt^ strategy [T7J (see the anti-correlation 
between Af ir and A s ± ix in Table 1 of [20]). 

The accuracy on these CP-asymmetries will be substantially improved 
once a precise measurement on A^ r will be available or the error of BR(Bd — > 
K°K°) and the QCD uncertainties on (mainly m c /m& and scale depen- 
dence) will be reduced. 



3.3. Supersymmetry 

The leading gluino-squark box and penguin contributions [3] to B s —* 
K°K° and B s — ► K + K~ were evaluated first in [18] using the U-spin flavour 
strategy with B^ — * 7r + 7r~ and, afterwards, in [27] using the new method 
combining flavour and QCDF [20]. The relative size of this contribution 
compared to the SM penguin is (a s /M^ usy )/(a/M^) ~ 1. The amplitude 
of these decays in presence of NP contains an extra contribution: A(B® — > 
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K+K-) = Af M + A u e^, A(B° S -» K°K°) = Af M + A d f**. These NP 
amplitudes A u e^ u and _4. rf e* d in terms of Wilson coefficients are: 



V2 



x(\4 + 4) - \(4 - 4) - (4 - 4) - ^cf g {i + 1; 



A 



with q = u,d, A = i(m 2 B - m 2 K )fKF Bs ^ K and x = 1-18 (see [27] for defi- 
nitions). These Wilson coefficients are sensitive to the s — b mass splitting. 
After including the constraints coming from BR(B — > Xg'j), B — > 7ri^ and 
AMg we found that a large isospin violation controlled by the mass splitting 
ur-cLr is possible between the NP amplitudes A u e % ® u and A d e l ® d . For the 
region of parameters considered in this super symmetric scenario, A u e i ^ u 
can be up to a factor three larger than A d e % ^ d . The specific results in 
supersymmetry for each decay mode are [27] : 

• Bg — > K + K~: The branching ratio is very little affected by SUSY. At 
most, the SM prediction can be increased by 15% for A^ = 0.1, in- 
creasing a bit the already good agreement with the new CDF data [8j . 
The direct CP asymmetry within SUSY falls inside the range —0.1 < 
A dir (B° s -> K + K~) SUSY < 0.7 for -0.1 < Af ir < 0.1. The deviation 
depends on the relative size of the competing SM tree versus the NP 
amplitude. A m i x {B® — > K + K~) SUSY can take any value from [-1,1]. 

• B° s -> K°K°: The impact of SUSY on BR(B Q S K°K°) is even 
smaller, reflecting the reduced allowed region for A d e i ^ d as compared 
to A u e^ u . The situation is very different for the CP asymmetries, 
that are particularly promising, due to the very small size of their SM 
prediction [20]. The direct CP asymmetry in SUSY can be 10 times 
larger than the SM one. A m i x {B® — ► K°K°) SUSY covers the entire 
range, and so this asymmetry can be large in the presence of SUSY, 
contrary to the SM prediction. 

The method discussed here is being applied to other non-leptonic -B-decays. 
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